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摘要 : 以 线粒体 细胞 色素 9 (Cyt b). 基因 作为 分 子 标 记 ， 对 晴 科 晴 亚 科 3 种 、 益 晴 亚 科 2 种 、 荔 晴 亚 科 1 
种 、 盾 晴 亚 科 2 种 昆虫 进行 序列 测定 ， 获 得 Cyt b 基因 432 bp 的 序列 片段 ， 该 片段 中 碱 基 T、C、A、6G 的 平均 
& 584 9 31.395. 12.496, 37.695580 18.7%，A + 平均 含量 为 68.9%， 明 显 高 于 G+C 含量 (30.190); 密 
码 子 第 三 位 点 A + TT 含量 更 高 达 82.7% 。 属 和 种 间 序 列 变异 大 ， 碱 基 替 换 多 发 生 在 第 三 位 点 。 以 得 豆 龟 晴 
( Megacopta cribraria) 为 外 群 构建 系统 发 育 树 ， 结 合 形态 特征 与 序列 变异 率 ， 先 同 将 盾 晴 亚 科 、 荔 晴 亚 科 从 晴 科 
划分 出 来 并 提升 为 科 的 观点 ， 它 们 之 间 的 系统 关系 为 盾 晴 科 与 荔 晴 科 形 成 姊妹 群 ， 较 晴 科 发 育 得 早 ; 晴 科 作 
为 一 个 单 系 群 ， 是 暑 总 科 中 最 为 进化 的 类 群 。 


xig: SES Cy 5b 基因; 系统 发 育 关 系 
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Discussion on the Phylogenetic Relationships of Partial 
Pentatomidae Insects Based on Sequences 
of Cytochrome 5 Gene 
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2. Life Science College, Shaanxi Normal University , Xi'an 710062, China) 


Abstract: This paper dealt with the nucleotide composition and the phylogenetic relationship among 8 species in 
Pentatomidae (Pentatominae, Asopinae, Tessaratominae, Scutellerinae) based on 432 bp sequences of Cyt b gene. The 
results showed: the average contents of T, C, A and G are 31.396, 12.496, 37.695 and 18.795, respectively; the 
contents of À T (68.946) is obviously higher than that of C + G (30.196), and come up to 82.746 at the third position 
of codon. The divergence of sequences is higher among genera and species, and replacement of nucleotides mostly occurs 
at the third position. The bug, Megacopta cribraria regarded as outgroup, the phylogenetic tree is set based on molecular 
information. Altogether with morphology characters, it is very reasonable that the two subfamily, Tessaratominae and 
Scutellerinae should be separated from family Pentatomidae, and be two valid families. The phylogenetic relationship of 
them is shown as follow: Tessaratomidae and Scutelleridae constitute a sister group and diverge more early than Penta- 
tonidae; as a monophyletic clade, Pentatomidae is the most evolutionary group in Pentatomoidea. 


Key words: Pentatomidae; Cyt b gene; Phylogenetic relationship 


KPREL (Pentatomidae) 是 半 翅 目 异 翅 亚 目 中 一 对 该 类 昆虫 的 系统 学 研究 ， 长 期 以 来 主要 集中 于 形 
个 较 大 的 类 群 ， 全 世界 已 知 760 属 4 100 种 左右 。 态 学 和 生态 学 方面 ,也 有 学 者 对 部 分 种 类 从 细胞 学 
(Schuh & Slater，1995)。 其 中 有 很 多 种 类 是 重要 的 ”等 方面 进行 了 研究 (Zhang & Zheng, 2001a, b; 
农林 业 害 虫 ， 严 重 危害 着 农林 业 和 畜牧 业 的 发 展 。 Ren, 2001; Zhao et al，1994)。 然 而 ， 由 于 这 一 类 
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群 种 类 多 ， 分 布 广 ， 且 许多 为 单 种 属 ， 再 加 上 形态 
构造 上 的 差异 ， 仅 依据 形态 性 状 对 其 进行 的 种 间 和 
亚 科 间 的 分 类 一 直 存 在 着 争议 。 国 外 许多 学 者 曾 将 
J*h. A58. Upéà 3 个 类 群 视 为 赌 科 下 设 的 亚 科 
(Miller, 1956; China & Miller, 1959), 我 国学 者 
Zhang (1985), Xiao (1977) 等 也 将 我 国 的 晴 科 昆 
虫 分 为 6 个 亚 科 ， 其 中 包括 上 述 3 个 亚 科 ， 且 这 一 
分 类 系统 在 国内 仍 被 沿用 。 近 年 来 ， 较 多 学 者 根据 
形态 构造 上 的 一 些 特点 ， 将 它们 从 晴 科 中 划 出 ， 分 
别提 升 为 科 级 水 平 (Schuh & Slater, 1995; Gapud, 
1991)。 同 时 将 晴 亚 科 中 的 一 些 属 划 出 ， 另 立 为 新 
的 亚 科 。 这 些 分 类 系统 均 是 基于 形态 学 内 容 建立 起 
来 的 ， 从 分 子 水 平 对 该 类 群 进行 系统 学 的 研究 国内 
鲜 见 报道 。 

线粒体 DNA (mitchondrial DNA, mtDNA) R 
有 严格 的 母系 遗传 ， 进 化 速率 较 核 DNA 快 ， 常 在 
一 个 种 的 存在 时 间 内 即 能 形成 可 用 于 系统 发 育 分 析 
的 性 状 状态 (Zhang & Shi, 1992; Huang，1998)。 
用 mtDNA 作为 分 子 标记 ， 已 广泛 用 于 昆虫 系统 学 
研究 中 (Chen & Zhang, 1999; Shao & Mao, 2002; 
Ren et al, 2002), fE mtDNA 中 ， 和 蛋白 质 编码 基因 
细胞 色素 b (cytochrome b, Cyt b) 基因 进化 速度 
适中 ,适合 种 间 甚 至 科 间 的 系统 发 育 研 究 (Yang et 
al，2001)， 是 目前 结构 和 功能 研究 得 最 为 清楚 的 
基因 之 一 (Gray，1989)。 本 研究 以 Cyt 基因 作为 


分 子 标 记 ， 首 次 对 晴 科 A 亚 科 8 种 昆虫 进行 特异 性 
扩 增 和 序列 测定 ， 并 以 龟 晴 科 的 筛 豆 龟 晴 (Mega- 
copta cribraria) 作为 外 群 ， 构 建 分 子 系统 树 ， 以 期 
从 中 获得 它们 的 系统 关系 ， 为 深入 研究 晴 科 昆虫 积 
累 分 子 生 物 学 的 证 据 和 资料 。 


1 材料 和 方法 


1.1 供 试 材料 

本 实验 所 用 昆虫 标本 均 为 2002 一 2003 年 采集 
的 新 鲜 标 本 ， 经 无 水 乙醇 固定 后 于 - 20 乞 保 存 。 
供 试 材料 种 类 、 采 集 日 期 、 地 点 见 表 1。 

1.2 DNA 的 提取 、 扩 增 和 测序 

剪 取 每 种 昆虫 一 个 个 体 前 胸 背 板 下 的 肌肉 少 
VF, EA 1.5 mL 离心 管内 ， 加 入 600 uL 匀 浆 液 
[A 液 :Tris 0.05 mol/L,NaCl 0.1 mol/L,Na; EDTA 
0.1 mol/L, pH 7.0- 8.0; B ii: 596 SDS; C 液 : 
人 蛋白酶 K 液 (2g/L); V (A) :V (B) :V (C) = 
8:1:1], 55 % 水 浴 3h, FHi/SUD hl RS, 10096 
的 冰 乙 醇 沉 淀 DNA, 7090 乙醇 洗 深 ， 室 温 干 燥 后 ， 
加 灭 菌 的 三 蒸 水 溶解 ，4 气 冰 箱 保 存 备用 。 

用 一 对 通用 引物 扩 增 线粒体 Cyt 5 基因 片段 ， 
引物 的 设计 参考 Simon et al (1994)， 引 物 序列 为 
CB1:5'-TATGTACTACCATGAGGACAAATATC-3' , 
CB2 : 5'- ATTACACCTCCTAATTITATTAGGAAT-3' , 

每 一 样品 的 扩 增 体积 为 50 pL, 反 应 体系 为 :10 x 








表 1 供 试 材料 名 称 及 来 源 
Table 1 Locations and classifications of tested samples 
亚 科 名 称 种 名 代号 采集 时 间 采集 地 点 数量 
Subfamily Species Code Collection time Locality Number 
1. *& VfL Pentatominae 菜 晴 Eurydema dominalus 
«T. | ,$ 
(Scopoli, 1763) Eul 2002. 7. 1l 四 川 万 源 Wanyuan, Sichuan 15 
横 纹 菜 晴 Eurydema gebleri 
E C i " 1 
Kolenati, 1856 u2 2002. 7. 1 陕西 黑龙 口 Heilongkou ，Shaanxi 9 
S8 Palomena angulosa pi j 
Motschulsky, 1861 Pal. a 2002. 9. 3 陕西 红河 谷 Honghegu, Shaanxi 17 
2 亚 科 Asopi Pinthaeus h li 
et Pim. h ”2002. 9. 11 。 陕西 杆 水 Zhashui，Shaanxi 7 
Horvath, 1911 
d$ Picromerus lewisi ; "OM 7 
Scot, 1874 Pir. 1 2003. 7. 16 ”陕西 平 利 Pingli，Shaanxi 7 
3. SH WB 硕 晴 Eurostus validus 9 ! ! 
Te dnm Dallas, 1851 Eu. v . 2002. 8. 26 ”陕西 南 五 台 Nanwutai, Shaanxi 3 
4. JE *É WE PE Scutellerinae — JE 8 Eurygaster testudinarius . ! 
(Geoffroy, 1785) Eyg. 2003. 7. 23 陕西 镇 坪 Zhenping, Shaanxi 8 
金 绿 宽 盾 晴 Poecilocoris lewisi . . 
(Distant, 1883) Poe. 2003. 7. 23 陕西 镇 坪 Zhenping, Shaanxi 6 
龟 晴 科 〈 外 群 ) 豆 龟 晴 M. ibrari 
第 TUE NOR outg 2003. 7. 24 四 川 万 源 Wanyuan, Sichuan 8 


Plataspidae (outgruop) (Fabricius, 1798) 
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buffer 5.0 uL, 25 mmol/L 的 MgCl; 4.0 uL, 0.25 
mmol/L dNTPs 1.25 pL, 10 pmol/uL 的 一 对 引物 各 
1.5uL, £j 50 ng/uL 的 模板 DNA 2.0 uL, TaqDNA 
聚合 酶 1.5U， 加 入 无 菌 三 蒸 水 至 50 忆 。 扩 增 条 件 
为 94 "C AE TE 4 min 后 ， 按 下 列 参 数 进行 35 个 循 
环 : 94 乞 变 性 30S, 50 乞 退 火 1min，72 延伸 1 
min， 最 后 一 个 循环 后 ，72 CBE $ min。 扩 增产 
物 用 1.2% 的 琼脂 糖 凝 胶 电泳 检测 ， 对 扩 增 效果 良 
好 的 样品 委托 上 海 生 工 生物 工程 有 限 公 司 进行 纯化 
和 测序 ， 测 序 仪 为 ABI PRISM377 型 。 为 保证 序列 
的 正确 性 ， 所 有 样品 均 进行 双向 测序 。 
1.3 DNA 序列 数据 处 理 

将 每 对 序列 用 ContigExpress 软件 进行 正 反 链 的 
拼接 检测 ， 并 与 序列 图 仔细 校对 ， 使 所 测序 列 保证 
为 两 特异 性 引物 之 间 的 原始 序列 。 用 Clustal X 
(1.8) 软件 对 序列 进行 比 对 。 用 分 子 进化 遗传 分 析 
软件 MEGA2.1 统计 序列 转换 、 苏 换 数 ， 分 析 各 物 
种 间 的 序列 差异 。 采 用 邻接 法 (Neighbor-Joining 
method，NJ )、 最 简约 法 ( Maximum Parsimony 
method, MP) 构建 系统 树 。 同 时 用 PHYLIP3.6a 软 
件 包 中 的 极 似 然 法 (Maximum Likelihood method, 
ML) 建树 。 


2 结 果 


2. Cyt b 基因 序列 组 成 及 变异 

本 实验 测定 了 晴 科 4 亚 科 8 种 及 外 群 共 9 条 线 
粒 体 DNA Cyt b 基因 中 间 的 部 分 序列 ， 长 度 为 433 
bp (不 包括 引物 长 度 ) ， 通 过 与 GeneBank 中 的 Cyt 





b 基因 全 序列 进行 对 位 排列 ， 得 知 所 测序 列 的 第 一 
个 碱 基 是 前 一 个 氨基 酸 密码 子 的 第 三 位 点 ， 为 分 析 
方便 ， 将 第 一 个 碱 基 去 除 ， 仅 以 432 bp 长 度 的 序列 
进行 分 析 ， 序 列 中 无 碱 基 的 插入 和 缺失 。8 条 序列 
中 (不 包括 外 群 ) 有 182 个 核 甘酸 变异 位 点 ，115 
个 简约 信息 位 点 。 碱 基 T、C、A、G 的 含量 分 别 为 
31.396, 12.496. 37.696 $0 18.726, A- TOES 
量 为 68.9% ， 明 显 高 于 G +C E (30.190), 
每 个 氨基 酸 密码 子 来 看 ， 碱 基 组 成 在 不 同位 点 差异 
很 大 ， 第 三 位 点 有 很 高 的 A+T 人 含量 达到 
82.7% 。 在 氨基 酸 组 成 上 ， 共 编码 144 个 氨基 酸 ， 
其 中 有 71 个 发 生变 蜡 ， 变 异 数 占 49.8%。 

统计 两 两 序列 间 的 碱 基 替 换 数 和 序列 差异 百 分 
比 〈( 表 2)。 总 体 上 看 ， 序 列 间 转 换 (transition) 略 
多 于 颠 换 (transversion) ， 且 碱 基 蔡 换 多 数 发 生 在 密 
码 子 第 三 位 点 ， 占 总 数 的 69.6% 。 
2.2 系统 树 的 构建 

以 Megacopta cribraria 作为 外 群 构 建 系统 发 育 
Bi, SET Kimura-2 参数 模型 ， 选 取 转 换血 换 信 
息 ， 同 时 应 用 自 举 检测 (bootstrap test) 以 估计 各 结 
点 的 自 引导 值 (bootstrap confidence level), 

由 系统 发 育 树 看 出 ， 属 内 种 间 的 个 体 首先 以 较 
高 的 置信 和 度 聚 在 一 起 ， 然 后 与 同一 亚 科 的 属相 聚 ， 
表明 种 间 差 异 小 于 属 间 差 异 小 于 亚 科 间 差 异 。 在 NJ 
树 (图 1A) 中 ， 晴 亚 科 的 2 属 形成 姊妹 群 后 先 与 
益 晴 亚 科 的 并 晴 相 聚 ， 之 后 与 益 晴 以 很 高 的 置信 值 
相聚 形成 一 支 ; 盾 晴 亚 科 的 2 属 形成 姊妹 群 ， 尽 管 
置信 值 不 高 ,但 表明 二 者 间 的 关系 较 其 他 属 近 ,这 


表 2 9 种 晨 象 碱 基 转 换 数 / 颠 换 数 (下 三 角 ) 和 序列 差异 百分比 (上 三 角 ) 
Table 2 The numbers of nucleotide transitions/transversions (below diagonal) and percents of sequences divergence 
(above diagonal) in nine species of pentatomid bugs 





1 2 3 4 5 6 7 8 9 

1 Eul 0.090 0.120 0.174 0.183 0.220 0.234 0.192 0.243 
2 Eu2 33/6 0.127 0.178 0.169 0.206 0.213 0.176 0.236 
3 Pal.a 24/28 21/28 0.160 0.181 0.218 0.222 0.204 0.241 
4 Pin.h 36/39 42/35 36/33 0.190 0.238 0.236 0.229 0.238 
5 Pir.1 48/31 48/25 41/37 50/32 0.220 0.225 0.218 0.241 
6 Poe.1 43/52 39/50 40/54 56/41 46/49 0.220 0.218 0.252 
7 Eyg.t 56/45 41/45 49/41 54/48 49/48 52/43 0.236 0.229 
8 Eur.v 43/40 38/38 42/46 50/49 41/41 42/52 51/51 0.220 
9 outg 56/49 55/41 55/49 51/46 52/52 60/49 51/42 54/41 


1 一 9 代表 种 类 同 表 1 (1 -9 present the same species as Table 1)。 
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图 1 9 种 晴 象 Cyt 5b 基因 序列 片段 NJ BERI MP 树 


Fig.1 NJ tree and MP tree based on Cyt 5 gene sequences fragment in nine species of pentatomid bugs 


代码 代表 种 类 同 表 1 (Code present the same species as Table 1)。 


一 分 支 与 第 一 支 接近 。 硕 晴 则 表现 为 单 系 性 ， 位 于 
系统 树 的 基部 ,是 分 化 较 早 的 一 支 。MP 树 (图 1B) 
中 ,部 分 属 种 的 位 置 发 生 了 变化 ,用 PHYLIP3.6a 
软件 包 中 ML 法 建树 ， 其 拓扑 结构 与 MP 树 完全 一 
致 。 
3 i it 

对 所 测序 列 分 析 的 结果 表明 ， 晴 科 8 种 昆虫 的 
A+ 人 TT 平均 含量 (68.9%) 明显 高 于 G+cC 含 量 
(30.1% )。 虽 然 比 绢 时 (75.4%)、 意 大 利 蜜蜂 
(84.995) TRURSà (77.490) 等 昆虫 线粒体 基因 组 
的 A+T 含 量 均 低 ， 但 仍 高 于 昆虫 A+T 的 平均 值 
(6496) (Chen & Zhang, 1999; Crozier & Crozier, 
1993; Jermiin & Crozier，1994)。 而 第 三 位 点 的 A 
+ 了 含量 高 达 82.7% ， 反 映 出 Cyt b 基因 在 密码 子 
的 使 用 上 具有 明显 的 偏向 性 。 从 碱 基 蔡 换 的 结果 
看 ， 密 码 子 第 三 位 点 的 替换 频率 显著 高 于 前 两 个 位 
点 ， 这 与 密码 子 第 三 位 点 承担 较 小 的 进化 压力 ， 发 
生 的 碱 基 蔡 换 较 少 引起 氨基 酸 的 替代 相 一 致 。 转 换 
的 发 生 主要 以 A<>G 为 主 ， 烟 换 的 发 生 主 要 以 A 
一 >T 为 主 ， 其 他 类 型 的 替换 很 少 发 生 。Crozier & 
Crozier (1992) 通过 对 高 A +T 与 氨基 酸 序列 关系 
的 广泛 研究 ， 认 为 高 A +T 与 氨基 酸 序列 变异 没有 
显著 的 相关 性 ; Desalle et al (1987) 则 认为 高 A+ 
T 增 加 了 AT 颠 换 的 可 能 性 ， 因 而 导致 第 三 位 点 的 
颠 换 率 增高 ; Hugall et al (1997) 也 认为 高 A+T 
加 速 了 氨基 酸 序列 的 变异 。 本 研究 所 揭示 的 晴 科 8 
种 昆虫 Cyt b. 基因 序列 中 的 高 A € T 与 氨基 酸 序列 
变异 之 间 的 关系 还 有 待 于 进一步 深入 地 研究 。 


Knight & Mibdell (1993) 认为 在 重建 系统 树 
时 ， 转 换 与 颠 换 比 的 值 若 小 于 2.0， 则 此 基因 序列 
突变 已 达到 饱和 状态 ， 受 进化 噪音 的 影响 可 能 性 较 
大 ， 重 建 系统 发 生 时 如 不 进行 特别 加 权 就 会 出 现 错 
误 信息 。 本 实验 所 获得 的 9 条 序列 转换 与 颠 换 比值 
平均 为 1.2， 为 估计 碱 基 替 换 的 饱和 程度 ， 以 遗传 
距离 (P-distance) 为 横 坐 标 ， 转 换 和 颐 换 为 纵 坐 标 
作 散 点 图 (图 2)。 并 依据 2 种 分 子 模型 Kimura-2 
参数 模型 、Jukes-Cantor 参数 模型 计算 种 间 遗 传 距 
离 ( 表 3)， 得 到 的 2 种 遗传 距离 矩阵 反映 的 遗传 距 
离 基本 一 致 ， 尽管 第 三 位 点 无 义 替换 高 于 第 一 、 二 
位 点 , 但 由 图 2 可 以 看 出 ， 转 换 和 颐 换 与 遗传 距离 
之 间 有 着 良好 的 线性 关系 ， 故 在 序列 比 对 过 程 中 对 
转换 、 颠 换 及 密码 子 位 点 未 作 加 权 ,而 全 部 采用 等 


70r Qo 转换 Transition 
C) 颠 换 Transversion 9 


转换 数 ( 颠 换 数 ) 


No. of transition (No. of transversion) 





0 0.1 0.2 0.3 
遗传 距离 P-distance 


图 2 9 种 晴 象 细胞 色素 b 基因 部 分 序列 替换 饱和 
度 分 析 散 点 图 
Fig.2 Variation of substitution versus P-distance of Cyt 
b sequences fragment of 9 species bugs 
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表 3 9S EE (下 三 角 为 Jakes-Cantor 距离 ， 上 三 角 为 Kimura 2 距离 ) 
Table 3 Pairwise distance of Jukes-Cantor-parameter (below diagonal) and Kimura 2-paramete (above diagonal) 


in nine species of pentatomid bugs 





1 2 3 4 5 6 7 8 9 
] Eu! 0.098 0.131 0.198 0.213 0.261 0.284 0.224 0.297 
2 Eu2 0.096 0.140 0.205 0.195 0.242 0.253 0.202 0.287 
3 Pal.a 0.131 0.140 0.181 0.208 0.258 0.266 0.239 0.294 
4 Pin.h 0.197 0.204 0.180 0.222 0.291 0.287 0.276 0.291 
5 Pir.l 0.210 0.191 0.207 0.219 0.262 0.269 0.259 0.293 
6 Poe.l 0.260 0.241 0.257 0.287 0.260 0.263 0.258 0.312 
7 Eyg.t 0.280 0.251 0.264 0.284 0.267 0.260 0.286 0.278 
8 Eur.v 0.222 0.200 0.238 0.273 0.257 0.257 0.284 0.264 
9 outg 0.294 0.284 0.290 0.287 0.290 0.308 0.273 0.260 
1 一 9 代表 种 类 同 表 1 (1 — 9 present the same species as Table 1)。 


权 处 理 。 

从 系统 树 看 ， 晴 亚 科 与 益 晴 亚 科 关系 最 为 接 
近 ， 而 且 是 晴 科 中 最 为 进化 的 类 群 ， 这 与 传统 的 分 
类 观点 一 致 。 根 据 序 列 变异 ， 拷 晴 亚 科 、 荔 晴 亚 科 
与 前 二 者 种 属 间 的 序列 差异 明显 要 高 ， 最 高 达到 
23.8% 和 22.9% 。 在 现行 的 分 类 系统 中 ， 多 数学 者 
则 意 将 蔓 晴 亚 科 、 盾 晴 亚 科 划 出 晴 科 另 立 为 科 : 从 
形态 特征 看 ， 盾 晴 前 胸 侧 板 向 前 扩展 成 游离 的 叶片 
状 ， 中 胸 小 盾 片 发 过 ， 盖 住 整个 腹部 和 前 起 的 绝 大 
部 分 ， 第 二 腹 节 气门 外 露 CREER ICE), qi 
晴 科 中 的 其 他 类 群 没有 上 述 特征 ; 此 外 ， 二 者 的 受 
精 赛 构 造 与 晴 科 中 其 他 昆虫 不 同 ， 晴 科 昆 虫 雌性 受 
TidtdE EET IPESEKXI BUDLBSERME , 
J& VA RED AGIT AL EEUE RODUE RI RE IRLRE EM 
段 骨 化 的 受精 囊 管 包围 其 中 ， 这 一 特征 是 晴 科 昆虫 
仅 有 的 ， 而 盾 晴 、 荔 晴 均 无 此 特征 (Zheng & Gui, 
1998)。 尺 管 目前 还 没有 资料 显示 晨 科 昆虫 中 序列 
差异 数据 与 分 类 单元 的 确立 之 间 存在 着 密切 的 相关 
性 ， 但 通过 本 实验 结果 ， 结 合 形态 特征 ， 我 们 亦 赞 
同 式 晴 亚 科 、 盾 晴 亚 科 从 晴 科 中 分 立 出 来 并 提升 为 
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